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The evolution of a large-scale azimuthal velocity field in a rotating cylindrical layer of fluid (radius 150 mm, depth 30 mm, free upper 
surface) with meridional convective circulation was studied experimentally. Two cases were considered: direct circulation provided by 
a rim heater at the periphery and indirect circulation provided by a central heater. The heating rate is characterized by the Grasshoff 
number Grf defined through the density of the heat flux. The detailed 3D structure of the mean large-scale velocity field is reconstructed 
using the PIV technique for 10 5 < Gr f < 4 - 10 7 . It was shown that the energy of meridional circulation grows with the Grasshoff number 
as i/Grf in both directions of circulation. Due to the action of the Coriolis force the meridional flow results in differential rotation. 
Differential rotation is characterized by the mean values of radial and vertical gradients of azimuthal velocity. Strong negative mean 
radial gradient which grows with the Grasshoff number is provided by direct circulation. In the case of indirect circulation a pronounced 
negative gradient arises at moderate Grashoff number. The behavior of the mean vertical gradient is quite different: a positive vertical 
gradient grows logarithmically with Grasshoff number under direct circulation, whereas a weak negative gradient characterizes the indirect 
circulation. This difference follows from the structure of the flow - the direct circulation provides a large cyclonic area localized above 
the anticyclonic flow, while indirect circulation leads to a strong separation of these two areas in the radial direction (the central part 
is occupied by the cyclonic flow and the periphery by the anticyclonic flow). Meridional circulation leads to substantial variation of the 
integral angular momentum. Direct circulation results in the growth of the integral angular momentum and indirect circulation causes 
it to decrease. At the same heating power, the increase of angular momentum at direct circulation is much stronger than its decrease at 
indirect circulation. 



1 Introduction 



All cosmic bodies rotate to a greater or lesser extent. The evolution of these bodies (excluding catastrophic 
events) proceeds under a very slow variation of the global angular momentum. However, if the cosmic 
object is gaseous or liquid (galaxies, stars) or includes liquid and/or gaseous shells (planets) the averaged 
large-scale azimuthal velocity field becomes essentially inhomogeneous in both the radial and meridional 
directions. The deviation from the solid body rotation is called diff erential rotation (PR) . 

DR plays a crucial role in the generation of cosmic magnetic fields ( Zeldovich et al. 19831 ) and has created 
considerable interest in investigating DR in electri cally conducting cosmic media - such as stellar convec- 
tive z ones, liquid planet cores and galactic disks ( Kleeorin &, Rogachevskii 20061 . Kitchatinov &: Rudigerl 
20071 ) . The study of DR in nonconducting rotating (spherical) layers has primarily been motivated by 
atmospheri c research beca use DR is a part of the global atmospheric circulation that determines climate 
formation ( Williams] 19681 ). For the most part, large-scale circulation in the atmosphere is caused by the 
horizontal temperature gradient, i.e., it has a convective nature . This strongly motivated laboratory exper- 
imen ts on convection in rotating vessels (cylinders or annuli) (jrlighnett et al\ 1981 , Boubnov &; Golitsvnl 
1995j). The first generation of laboratory experiments primarily intended to reproduce the Hadley cell (or 
direct meridional circulation) by realiz ing a conv e ctive flow in rotating annuli or flat cylinders heated at 
the periphery and cooled in the center (Hidelll958l . Fultz Kaiser 197ll ). Wave-transition spectra and vac- 
illations were studied bv lSpence fc Fultzj (|l977n T Interest in DR in this configuration has b een heightened 
by the discovery of so-called super-rotation in the atmosphere of Venus ( Belton et al. 19761 ) - the existence 
in the high atmosphere of very fast latitudinal circulation, which exceeds by about 60 times the rotational 
velocity of the surface (the rotation period is about four days versus 243 days of the Venus astronomi- 
cal day). The problem of super-rotation is directly linked with processes of angular momentum diffusion 
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and transport in planetary atmospheres. It was proposed by Gierasch ( 19751 ) that Venus super-rotation is 
maintained by angular momentum transport provided by meridional circulation. S ome mechanisms, which 
may contribu te to the development of a strong differential rotation were proposed ( Shubert Young 1970 . 
Rossowill983r). Compreh ensive review of different quasy-axisimmetric models of planetary atmospheres can 
be found in (IR,eadlll986h . 

The first experimental attempt to quantitatively analyze meridional transport of ang ular momentum 
i n the rotating layer of fluid in the case of direct meridional circulation was realized by ISpence &: Fultz 
however, integral characteristics such as overall angular momentum were not studied. Numerical 
simulations for a similarly stated problem using different boundary condit ions were presen t ed by IWilliamd 
(1967, 19681 ) and showed good agreement with experimental results from lSpence Fultd (|1977V ). 

Properties of angular momentum budget and super-rotation of axis ymmet r ic the rmally-driven circula- 
tions in a rotating cylidrical fluid annulus were studied numerically by Read ( 19861 ). Angular momentum 
diffusion due to molecular viscosity was examined and was found to be very important for super-rotating 
flow formation in a system with stress-free top and side boundaries and a non-slip bottom. It was shown 
that global super-rotation is connected with direct meridional circulation. 

An opposite situation to the direct meridional cell is the indirect cell, usually produced in a rotating 
layer of fluid by heating at the center. Mainly, it was studied in the context of another atmospheric 
problem concerning th e generation of intense large-scale vortices. Laboratory modeling of these typhoon- 
like vortices was made ( Bogatvrev 199Cll . Bogatvrev fe Smorodinlll996 . Bogatvrev et aljfeood ). These studies 
were concentrated on cyclonic vortex formation and its evolution. A qualitative experimental study of 
convective flow driven by a finite-sized circular heating plate at the bottom of a horizon tal fluid layer, 
both with and without background rotation, was carried out by Boubnov Heijstl ( 19941 ). The angular 
momentum budget for an indirect cell was analyzed by Read ( 19861 ) and it was shown that indirect 
circulation resulted in local and/or global sub-rotation, thereby demonstrating that the change in overall 
angular momentum essentially depends on the type of meridional circulation. 

In this paper, we return to experimental investigation of the convection flow in rotating cylindrical 
vessels, being motivated by following reasons. First, we examine the differential rotation problem in the 
rotating layer heated at the periphery or center from a uniform viewpoint, paying special attention to 
the integral characteristics of differential rotation. Second, our experiment concerns the shallow layer of 
fluid. And third, we obtain the full structure of mer i dional and azim uthal velocity fields exploiting the 



potentialities of the PIV technique ( Raffel et al. 1983 . Xia et al. 19761 ). 



We note that our statement of the problem is similar to that made in riR,eadl[l9 86). except for the aspect 
ratio and the physical properties of the fluid. Numerical simulations in (jReadlll986l ) were carried out for 
fluid layers in an annular channel with an aspect ratio close to unity, whereas we studied the cylinder layer 
(no inner wall) with a fixed aspect ratio e = h/R = 0-2 (h - depth of the layer, R - radius of the cylindrical 
vessel). The Prandtl number for the fluid in (|Readl [l986) was about 10, and in our study it is about 100. 

The structure of the paper is as follows. The experimental set up is described in Sec. El The structure 
of convective flows for direct and indirect meridional cells is represented in Sec. [3] and Sec. 01 respectively. 
Integral characteristics of DR are described in Sec. [5j Our results are summarized in Sec. [6) 



2 Experimental set-up 

We studied a convective flow in a flat cylindrical vessel placed on a rotating horizontal table(Fig.[T]). The 
table provides uniform rotation in the angular speed range 0.04 < Q < 0.30 s" 1 (with accuracy of ±0.002 
s _1 ). The cylinder radius is R = 150 mm. The cylinder is made of acrylic plastic, which is transparent to 
the laser sheet. The fluid is transformer oil, characterized by a high Prandtl number 

Pr = v/x, (1) 

where v, x are the kinematic viscosity and thermal conductivity (for transformer oil, Pr = 116 at 50° C). 
The layer thickness is constant in all experiments {h = 30 mm) and the upper surface is free. 
Two heaters were installed at the bottom of the vessel (Fig. [I]). Heater I is a copper rim 20 mm wide 



June 7, 2008 4:27 Geophysical and Astrophysical Fluid Dynamics DIF'ROT'GAFD 

Laboratory study of differential rotation in a convective rotating layer 3 



CCD Camera 




Figure 1. Left panel - scheme of the laboratory set-up; Right panel (view from above) - zone of horizontal velocity measurements 
where a typical (here inversed) PIV image is shown. O - axis of rotation. 




Figure 2. Sketch of meridional circulation, initiated by peripheral heater I (direct circulation - a) and by central heater II (indirect 

circulation - b). The heaters are shown by black boxes. 



Tabic 1. Experimental parameters 
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placed at the periphery of the cylinder. This heater provided the direct meridional circulation as shown in 
Fig.[2|a). Heater II is a copper cylinder of radius 52 mm, installed in the center of the vessel and coinciding 
with the axes of rotation. The heating near the rotation axes (polar heating) provides circulation in the 
opposite direction (indirect circulation, Fig.^b)). Both heaters operate from a DC power supply. The 
room temperature is kept constant, and cooling is provided by the heat exchange on the free surface. It 
takes about two hours to obtain a steady-state temperature regime. 

The velocity field measurements were made for steady state regimes by a 2D particle image velocimetry 
(PIV) system, "Polis", manufactured by the Institute of Thermophysics (Novosibirsk). A CCD camera is 
placed above the rotating vessel in the laboratory coordinate system as shown in Fig.[TJ The measurements 
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Figure 3. Grasshoff number versus heating power: squares — heater I, circles - heater II, C = 0.069s 1 . 

are carried out in a confined area of the vessel, but the rotation provided averaging for the entire hori- 
zontal cross-section. The area of velocity measurements is shown in Fig.[TJ PIV measurements are done 
at horizontal cross-sections for 4 < z < 28, mm in 2 mm steps (the cylindrical coordinate system (r, (ft, z) 
is used). We focused on flow motions in the rotating frame, and for this purpose the velocity field of 
solid-body rotation was subtracted and all velocity fields and flow characteristics, described as follows, are 
related to the motions in the rotating frame. The time-averaged velocity fields consist of 41x30 vectors 
and were averaged over 100 instant velocity fields. Then time-averaged horizontal velocity components 
and v r are averaged along the coordinate <j) and the mean velocity fields at vertical cross-section are 
reconstructed. The domain of PIV measurements is restricted by the boundary layers, where the velocity 
field reconstruction is hampered by optical distortions caused by strong temperature gradients near the 
heaters and/or by high velocity gradients. These boundary layers are shown in white in figures given below. 
The PIV velocity measurements were accurate to within 5%. 

Two main factors are responsible for DR in the flows under discussion - convection and rotation. To 
characterize the intensity of convective flows, we use the Grasshoff number Grf , defined in terms of layer 
depth h and heat flux density q = Pj S n (P is the power of the heater and S n is the heater's surface area) 



Gr f 



cpxv 2 ' 



(2) 



where g is the gravitational acceleration, is the coefficient of thermal expansion, c is the thermal capacity, 
and p is the density. Because the fluid characteristics depend on temperature (particularly viscosity) and 
the mean temperature of the layer is a function of the heating power P, Grf is a nonlinear function of P. 
In Fig. [3] Grf versus P is plotted on a log-log scale for both heaters. This graph shows that the dependence 
is even stronger than a squared relationship, and is similar for both heaters. 
As a non-dimensional characteristic of rotation we use the Ekman number 



E 



2Qh 2 



(3) 



because formation of DR strongly depends on angular momentum exchange in viscous boundary layers. 
The values of parameters for all experiments are given in Table[TJ 



3 Direct circulation 



In laboratory models of general atmospheric motions, the direct meri dional circulation is usually pro- 
vided by heating the outer (peripheral) p arts of cylindrical vessels ( Highnett et al. 1981 . Hide! 1958 . 
Fultz k Kaisei]ll97ll . ISoence k Fultzl Il977l ). The generated flow can be considered as a rough model of 
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Figure 4. (Color online.) Stream function of the meridional mean velocity field (a) and mean azimuthal velocity field (b) for direct 
circulation at Grf = 1.4 ■ 10 s , E = 0.118. Positive values of v$ are related to the cyclonic flow and negative- to the anticyclonic flow. 
The border between the cyclonic and anticyclonic motion is shown by the white line. 




atmosp heric meridiona l circulation (the Hadley cell) in which the role of latitudinal variation of VL Z is ne- 
glected ( Williams] 19681 ). In our experiments, a horizontal temperature difference is maintained by heating 



the periphery of the bottom and by natural cooling at the free upper surface. We studied the evolution 
of the direct meridional flow structure with variation in the heating power and the contribution of the 
meridional toroidal cell to differential rotation. All measurements are made for a fixed rotation speed, 
n = 0.069s- 1 . 

At weak heating the meridional flow has the structure shown in Fig.[2ja). The stream function ip, defined 
for the mean flow in the radial cross-section as 

d z tp = -rv r , d r i\) = rv z , (4) 

where v r is the radial velocity component and v z is the vertical velocity component, is shown in Fig.UJa for 
Gr f = 1.4 • 10 5 , E = 0.118. The meridional cell occupies the whole layer providing an inward (polarward) 
radial flow in the upper part and an outward flow in the lower part of the layer. The center of the cell 
(the maximum of the stream function) is localized close to the cylinder wall, where the heater produces 
an intensive upward flow. 

Formation of differential rotation in such a system can be described as follows. At first, the action of 
the Coriolis force on the radial flows initiates cyclonic (prograde) azimuthal flow in the upper layer and 
anticyclonic (retrograde) flow near the bottom. Then, angular momentum transport provided by meridional 
circulation and angular momentum diffusion due to molecular viscosity lead to the steady-state regime 
shown in Fig.[Ujb). The maximum cyclonic velocity is located near the free surface at r ~ R/2, and the 
maximum anticyclonic velocity is located near the bottom, shifted to the outer side. 

The growth of Grf intensifies the meridional flow (see Fig.[5ja), where the stream function for Grf = 
1.4 • 10 7 is shown). Now, the main meridional convective cell is pressed more to the periphery of the 
rotating layer. The azimuthal velocity field in Fig.[5](b) shows an absolute dominance of the cyclonic zonal 
circulation, and the anticyclonic motion survives only in a small area near the wall. The maximum of the 
cyclonic flow is shifted polarward compared to its position in the previous case (Fig.|4|), but the intensity 
of this flow grows dramatically (the maximal azimuthal velocity grows from 0.7mm/s to 13 mm/s). The 
ratio of the maximum of the cyclonic flow velocity to the velocity of the underlying solid border is more 
than 3 (compared to 0.1 in the previous case). 

The measured velocity fields are in good qualitative agreement with the numerical simulation of rotating 
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Figure 6. (Color online.) Stream function of the meridional mean velocity field (a) and mean azimuthal velocity field (b) for indirect 
circulation at Grf = 4.1 ■ 10 5 , E = 0.104. The dashed black line marks the border between the cyclonic and anticyclonic motion. 
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Figure 7. (Color online.) Stream function of the the meridional mean velocity field (a) and mean azimuthal velocity field (b) for 
indirect circulation at Grf = 3.02 • 10 7 , E = 0.032. The solid white line marks the border between the cyclonic and anticyclonic motion. 



convection in a wide annular gap ( Williams! 1968I ). The dominance of cyclonic flow in the layer with a free 
surface is mainly caused by different boundary conditions on the free surface and at the rigid bottom. I n 
the case with all rigid boundaries cyclonic and anticyclonic flow were essentially balanced ( Williams] 1968I ). 



4 Indirect circulation 

In this section we describe the mean flow induced by heater II, located at the center of the cylinder 
bottom. Natural convection over the local heat sources with and without background rotation has many 
geophysical and engineering applications and has been widely studied. Most of these studies were focused 
on heat transfer measurements and qualitative studies of the flow structure. A systematic study of the 
dependence of the c onnec tive flow regimes on the aspect ratio 5 and the Rayleigh number was carried out by 



Boubnov &: Heiistl (|1994l ). The aspect ratio was defined as 6 = D/h, where D is the diameter of the heater. 
Four domains corresponding to different regimes were defined in the (5, Ra) plan: I - laminar toroidal 
cell; II - thermal plumes; III - turbulent flow; and IV - transition regime. Following this classification, 
all our experiments belong to case III - turbulent flow (6 = 3.5, 1.5 • 10 5 < Ra^ < 6.5 • 10 6 ). Here 
Ra h = gPATh^ux)- 1 . 

The structure of the mean (axisymmetric) flow is determined again by two factors: heating power and 
frame rotation speed. First, we varied the heating power at a fixed rotational speed (Q = 0.069 s _1 ). Local 
heating at the center of the bottom generates vertical and horizontal temperature gradients. The horizontal 
temperature difference provides an indirect axisymmetric meridional cell (Fig.Efb)). The flow in the lower 
part is directed inward, tending to the central bottom region, where the heat source is located. A strong 
upward flux is formed at the center above the heat source. In the upper layer, the outward flow is directed 
toward the periphery. The fluid is cooled at the free surface and finally moves downward along the side 
wall. Fig. [6] shows the mean velocity fields at Grf = 4 • 10 5 . The growth of Grf increases the intensity of the 
meridional flow, but does not appreciably change its structure (compare panels (a) in Fig. [6] and Fig. [7]). 

Following the main direction of the meridional flow the cyclonic motion is accumulated near the rotation 
axis. At low heating the cyclonic flow occupies a small domain close to the bottom, and the anticyclonic 
flows dominates in the upper layer (see Fig.[6^b)). The increase of Grf leads to a substantial growth of 
the cyclonic velocity maximum (from 0.2 mm/s to 7 mm/s) and weak changes of anticyclonic velocity 
maximum (about 1 mm/s), compare panels (b) of Figs.[6l|71 The cyclonic flow at high values of Grf 
completely occupies the central part of the layer, spreading narrow wings along the bottom up to the 



June 7, 2008 



4:27 



Geophysical and Astrophysical Fluid Dynamics DIF'ROT'GAFD 

Laboratory study of differential rotation in a convective rotating layer 



7 




periphery. The anticyclonic flow is pushed towards the side wall. 

The second important factor defining the flow structure is the speed of the frame rotation. Its influence 
is illustrated in Fig. [HI where the velocity fields for the central part of the layer, mostly occupied by 
cyclonic vortex, are shown for Grf = 9.1 • 10 6 and two different values of the rotation speed - O = 0.046 
(E = 0.066) and O = 0.161 s~ 1 (E = 0.019). The growth of the frame rotation speed leads to some 
decrease in intensity of the meridional circulation (see panels (a) and (c) of Fig.[8j because the increase of 
the background rotation suppresses the vertical motions. The azimuthal velocity fields for these two O are 
shown in Fig.[8j(b,d). At a higher value of O the cyclonic flow is confined in the lower part of the layer. 



5 Integral characteristics of differential rotation 

The preceding results show that meridional circulation in both directions leads to a pronounced DR. The 
direct circulation initiates a strong cyclonic flow near the free surface at middle radii. The anticyclonic 
flow in this regime is localized near the bottom of the vessel and is weakened with increasing heating 
power. The result of indirect circulation is not quite symmetric. In this case, the cyclonic and anticyclonic 
flows are better separated in radial direction - the cyclone dominates in the central (polar) part of the 
layer, spreading at high values of heating power throughout the entire layer depth, while the anticyclone 
occupies the upper part of the periphery. 

To quantitatively estimate the balance between the cyclonic and anticyclonic motions, we calculated 
the integral angular momentum for the whole layer; more precisely, for the whole area available for PIV 
measurements. Specifically, we calculated the relative variation S of the integral mom entum L rel ative to 
the solid body rotation momentum of the layer L s defined as "global super- rotation" (|Readl [l986) 



L-L s 



where 



s = —r^, (5) 



ph—Si pR—84 r2ir 

L = dz rdr / rv^dcj), (6) 

JS t Js 3 Jo 

h—S 2 fR—Si 

3, 



L s = 2tt dz nr d dr, (7) 

J Si J S3 

and Si are the thickness of the boundary layers which are not accessible for PIV measurements. The 
results obtained for different values of Grf and a given rotation speed are shown in Fig.[9j The sign of the 
global angular momentum is directly related to the heating mode: peripheral heating provides the direct 
circulation that leads to the growth of the integral angular momentum of the fluid layer. On average, the 
fluid layer rotates faster than the background; therefore, this is a case of super-rotation. The maximum 
value of S achieved in the experiment is about 0.4. In contrast, central heating provides the indirect 
circulation that reduces the integral angular momentum; on average the layer rotates more slowly than 
the background (sub-rotation). Note that the increase in angular momentum at direct circulation is much 



June 7, 2008 4:27 Geophysical and Astrophysical Fluid Dynamics DIF'ROT'GAFD 

8 VLADIMIR BATALOV, ANDREY SUKHANOVSKY and PETER FRICK 

0.5 
0.4 
0.3 
0.2 
0.1 


< 

-0.1 
-0.2 



CO 



10 



10 



10 



10 



Gr, 



Figure 9. Global super-rotation S versus Grf for direct cell (squares) and indirect cell (circles) circulation, f2 = 0.069s 1 . 




Figure 10. (Color online.) Local super-rotation s: (a)— direct circulation for Grf = 1.4 ■ 10 7 , E = 0.032; (b)— indirect circulation for 

Gr f = 3.02 ■ 10 7 , E = 0.032. 



greater than its decrease at indirect circulation for the same value of the heating power. The lowest value 

of S in the case of the indirect cell was S ~ —0.16. 

Useful information might give the so-called "local super-rotation" s (|Readl[l~986l ). defined as 



v<f>r 
SIR 2 



1. 



(8) 



It shows the excess of angular momentum in comparison with highest possible angular momentum of 
the fluid element in a solid-body rotation state. The distribution of s for different types of meridional 
circulation is shown in Fig. [TO) This figure shows that even for high values of Grasshoff number the local 
super-rotation s is negative everywhere (sub-rotation). Note, that nume rical simulation in case of stress- 
free sidewalls showed that there is noticeable area withs > (|E,eadlll98fil b This proves the strong influence 
of the side wall boundary conditions on the distribution of angular momentum. 

Angular momentum redistribution in the fluid layer occurs due to the action of net torque due to friction 
at the bottom and side walls (the upper surface can be considered stress-free, and the torque on the upper 
surface is negligible). In the case of direct meridional circulation the anticyclonic flow in the lower layer 
generates torque that injects angular momentum into the layer. Meridional flow provides the transfer of 
angular momentum from the lower to upper part of the layer. The sink of the angular momentum occurs 
in the viscous boundary layers occupied by the cyclonic flow. In the case of an indirect cell the sink of 
the angular momentum is located in the central part of the bottom occupied by an intensive cyclonic 
vortex, and the source of the angular momentum is located at the upper part of the side wall - the area of 
anticyclonic flow. 

Thus, the flow structure and the imbalance between cyclonic and anticyclonic motion observed in the 
experiment is a consequence of different boundary conditions: non-slip on rigid bound aries, and stress -free 
on the surfa ce. This agrees with general conclusions made from numerical simulations ( Williams] 19681 ) and 
(|Readlll986h . 
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Tabic 2. Grasshoff number and the parameter Q for all experiments 



Gr f 


1.4- 10 & 


7.7- W b 


2.9 • 10 b 


5.8 • 10 b 


1.4 • 10 Y 
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0.0082 


0.0066 


0.0066 


0.0063 


0.007 



Note that some di screpances e xist between our experimental results and the analysis of the azimuthal 
circulation made in ( Read 1986), where the flow structure analysis was made through dimensionless pa- 
rameter Q which was defined using the aspect ratio, Ekman number and Rayleigh number Ra, 



Q = Ra~ 1/2 E- 



(9) 



Ra 



gf3L 3 AT 



(10) 



where L is the horizontal length scale and AT is the horizontal temperature drop. Q describes the ratio 
between thickness of thermal boundary layer (non-rotating) and that of the Ekman boundary layer. It was 
shown that global super-rotation depends on Q in the follow way: the super-rotation is almost constant 
for slow rotation (Q < (e/Pr) 2 ); it increases at moderate rotation ((e/Pr) 2 Q < 1), and it decreases quickly 
at rapid rotation (Q > 1). Calculating Q for all measured experimental regimes (taking L = R), we found 
that Q changes very weakly (see Table[2]). Thus, the whole variety of observed regimes exist under similar 
values of Q and belong to the case of "moderate" rotation. 

The independence of Q on Grf can be easily understood if we expand the expression for Q and rewrite 
it in the following form 



Q = 2Q 



h 



gPAT 



1 

pT 



1/2 



:n) 



In our experiments, all parameters from (jlip except Pr and AT were constant. The growth of Grf leads 
to an increase of the mean temperature of the fluid T and AT, but to the decrease in the value of Pr. The 
decrease of Pr with growth of T is caused by strong temperature dependence of oil viscosity v. As a resul t 
(PrAT) -1 / 2 and parameter Q are almost independent of Grf. Note that the simulations of Read (1986) 
were performed for e = 1, Pr = 10, a fixed temperature drop (heating) and varying Ekman number. In our 
experiments, the aspect ratio was smaller (e = 0.2), the Prandtl number an order of magnitude larger, the 
rotation speed fixed. The heating power was the main governing parameter, which provided the variety of 
regimes. 

As previously stated, the structure of the meridional flow is different at direct and indirect circulations 
(compare Fig.[5ja and Fig. 0a). To estimate the energy balance, we compared the total energy of the 
meridional flow W r for different kinds of circulation. Our measurements showed that the energy of the 
indirect meridional flow (provided by heater II) is about two times higher than the energy of the direct 
circulation (provided by heater I) for the same heating power. Therefore, for comparison of the energy 
budgets of the direct and indirect cells we used the Grasshoff number Grf and plotted the energy of the 
meridional circulation versus Grf. Fig. 1111 shows that despite different structures of meridional circulation 
for the direct and indirect cells W r for both cells is well correlated with the Grasshoff number. The function 
W r (Grf) is given in log-log scale and displays a power law close to "1/2". 

We then analyzed how the energy of the cyclonic and anticyclonic motions W c and W a depend on 
the type of circulation and on the Grasshof number. Direct circulation results in a much more intensive 
cyclonic flow - than indirect circulation at the same Grf. However, in spite of the great difference in the 
structures and intensities of the cyclonic flow in the direct and indirect cells the dependence of W c on the 
Grasshoff number is similar in both cases (see Fig.1121 where the energy of the cyclonic motion is shown 
versus the Grasshoff number for both circulations). 

We could not obtain a similar picture for anticyclonic flow, because in the case of direct circulation the 
anticiclonic motion is mainly located in the boundary layer and could not be reconstructed with sufficient 
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Figure 11. Total energy of the meridional flow versus Grasshoff number for direct (squares) and indirect (circles) circulation, 

f! = 0.069s -1 . The solid line shows the slope for W r ~ \/Gr. 
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Figure 12. Total energy of the cyclonic flow at direct (squares) and indirect (circles) circulation versus the Grasshoff number, 

Q = 0.069s" 1 . 



accuracy. Fig. [13] shows the ratio of the cyclonic to anticyclonic flow energy for the indirect circulation. 
From this graph it is evident that in the considered range of parameters, this ratio grows linearly with the 
Grasshoff number. Thus, if this tendency remains valid at higher values of Grj, the energy of the cyclonic 
flows will surpass the energy of the antiyclonic flow at Grj > 2 • 10 7 . 

Finally, we considered the direct characteristics of differential rotation. Note that DR is one of the 
major mechanisms of magnetic field generation in cosmic bodies; namely, DR is responsible for generation 
of the toroidal (azimuthal) magnetic field from the poloidal field. The efficiency of this mechanism depends 
on the gradient of the azimuthal velocity. The nondimensional characteristic of a dynamo process based 
on DR - the so-called "dynamo number" includes some characteristic value of the gradient of the mean 
azimuthal velocity (in the radi al or meridional direction) in the domain where the dynamo action is 



localized (jZeldovich et aZ.I Il983l ) . Since we studied the DR outside a special dynamo problem, we use, 



as characteristics of DR, the mean values of the gradients d r v$ and d z v^, calculated over the domain of 
reconstruction. 



rh — & 2 rR—&4 

D r = 2tt dz r(d r v<f,)dr, 

rh—S 2 rR—&4 

D z = 2ir / dz r{d z v ( f ) )dr. 
J Si J 5$ 



(12) 
(13) 



The results are shown in Fig.[T3]as functions of the Grasshoff number for both circulations. Panel (a) of 
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Figure 13. Ratio of the cyclonic motion energy to the energy of anticyclonic motion under indirect circulation, f! = 0.069s 1 . 
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Figure 14. Mean radial D r (a) and vertical D z (b) gradients of azimuthal velocity for direct (squares) and indirect (circles) circulation 

versus the Grasshoff number for, SI = 0.069s -1 . 



this figure shows that the mean radial gradient is negative for direct circulation (the central part rotates 
faster than the periphery), and that its value grows with the Grasshoff number. For indirect circulation, 
the mean radial gradient is close to zero at low Grasshoff number (the cyclonic and anticyclonic motions 
are separated vertically, but not radially) and becomes negative at high Grf . Then, the radial DR is weaker 
as in case of direct circulation, but the central part also rotates faster than the periphery. The behavior of 
the mean vertical gradient is quite different for two circulations: a positive vertical gradient increases with 
a growing Grasshoff number (approximately as log(Grf)) at direct circulation, whereas a weak negative 
gradient characterizes the indirect circulation. 



6 Summary 

We studied the formation of DR for direct and indirect meridional cells. In the case of a direct cell, the 
cyclonic flow prevails and spans nearly the entire layer except for a relatively small zone at the periphery 
near the bottom. The indirect cell is characterized by competition between the cyclonic flow in the central 
part and anticyclonic flow at the periphery. The ratio W c /W a grows with Grf and approaches unity at 
Gr f ss 2 • 10 7 . 

The structure of azimuthal flows and their intensities are determined by the difference between friction 
force torques at the rigid boundaries and the free surface. The condition of a steady-state regime is zero net 
torque. Torque due to frictional forces at the solid boundaries serves as a source of the angular momentum 
in the anticyclonic flow area and as a sink in the cyclonic flow area. During transition to a steady-state 
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regime the net torque is nonzero and leads to variation in overall angular momentum. Direct circulation 
leads to the growth of overall angular momentum, whereas indirect circulation causes it to decrease. At 
the same Grf the increase in angular momentum at direct circulation is much stronger than its decrease 
at indirect circulation. 

DR is caused by the meridional transport of angular momentum. We have shown that the energy of 
the meridional circulation grows with the Grasshoff number as \/Grf at both directions of circulation. We 
characterize the differential rotation by the mean values of radial and vertical gradients of the azimuthal 
velocity. Direct circulation provides a strong negative mean radial gradient, which means that the central 
part rotates faster. The same tendency exists at indirect circulation, but the gradient is much smaller or 
even close to zero at low Grasshoff numbers. Direct circulation gives rise to a pronounced vertical gradient 
(positive), which grows with Grasshoff number. Indirect circulation provides a weak negative gradient. 
This difference follows from the very structure of the flow: direct circulation provides a large cyclonic area 
localized above the anticyclonic flow, and indirect circulation leads to substantial separation of these two 
area in the radial direction (the cyclonic flow is localized close to the rotation axis and the anticyclonic 
flow is concentrated near the side wall). 
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